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We present a comprehensive study of the effect of a transverse magnetic field on the macroscopic 
quantum tunneling of two identical coupled particles with dissipation. Two coupled particles of 
identical masses but opposite charges in the plane, in the presence of a constant transverse external 
magnetic field and an external potential interacting with a bath of harmonic oscillators are studied. 
In this case, the problem cannot be mapped to a one-dimensional problem, it strictly remains two- 
dimensional. We calculate the exact effective action both for the case of linear coupling to the 
bath and without a linear coupling to the bath using imaginary time path integral and Leggett's 
prescription at T = 0. In the limit of zero magnetic field we recover a two dimensional version of the 
result found in Q for the case of two identical particles. We find that in the limit of weak and strong 
dissipation, the effective action reduces to a two dimensional version of the Caldeira-Leggett form 
in terms of the reduced mass and magnetic field. The case of Ohmic dissipation with the motion of 
the two particles damped by the Ohmic frictional constant n is studied in detail. 
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INTRODUCTION 

Macroscopic quantum tunneling with dissipation has 
become the subject of interest in quantum statistical me- 
chanics and condensed matter physics for many years 
EL 0, HI- This mainly involves the influence of the en- 
vironment (thermal bath of harmonic oscillators) on the 
tunneling of a macroscopic particle with variable say q 
out of an external potential V(q) which is assumed to 
have a metastable minimum. In most cases of physical 
interest, it is assumed that q interacts linearly with the 
environmental coordinate say x a (a = 1, 2, • • • ) at a cer- 
tain temperature T. The breakthrough in this subject 
was made by Caldeira and Leggett Q. They considered 
a Euclidean Lagrangian of the form 



C E = 



Mq 2 +V(q) + J^-m a (x 2 a 



2 2 \ 
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where the effective action is given by 
dr 



ceff 

Se = 





4tt 



hlq 2 + V(q) 
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and rj is the frictional constant. Chudnovsky 0| gen- 
eralized this formalism by considering two macroscopic 
particles that interact with each other via a nonlinear 
potential V(|a;i — 2:2!) with the coordinate x 2 linearly 
coupled to the environment. The Euclidean Lagrangian 
is of the form 



C E = \m x x\ + \m 2 x\ + V{\xi 
+ \-^2m a u 2 a (x a -x 2 ) 2 . 



X2\ 



\^m a x 2 a 



(G) 



where the parameters m a ,uj a ,c a need not to be known 
in detail. The partition function is given by 

K{q,x a -T) = j Vq(r) J Y[Vx a {r) cxp (S E ) , (2) 



where 



Si 



drCr 



(3) 



Performing the functional integral over x a in the limit 

t — ?• 00 gives 



K(q;r) = J Vq(r) exp (-5; 



E f 



(4) 



Integrating out the environmental degree of freedom and 
using the new coordinates 



xi - x 2 , 
M lXl + M 2 x 2 
Mi + M 2 



(7) 



he noticed that in the limit M\ —5- 00, the effective action 
reduces to the form of Caldeira and Leggett Eq.([5])- 

In this paper, we shall generalize Chudnovsky idea 
by considering two coupled macroscopic particles in the 
plane, in the presence of a constant transverse external 
magnetic field and an external potential. 
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MODEL 



We shall first consider the Euclidean Lagrangian 



c e = — Fi 

1 



2^2 



1 



Xi X X\ — X 2 X X 2 



-m 2 ulx2 2 + V(\xi - x 2 \ 



(8) 



This Lagrangian describes motion of two coupled parti- 
cles in the plane, in the presence of a constant transverse 
external magnetic field and an external potential. These 
two particles interact with each other by a nonlinear po- 
tential Vflicl — x 2 \) which has a metastable minimum. 
The magnetic field term in imaginary time corresponds 
to a symmetric gauge for the two particles. It breaks the 
time reversal symmetry of the Lagrangian and the weak 
harmonic oscillator potentials break the spatial transla- 
tion symmetry of the Lagrangian. however, we can re- 
store this symmetry (up to a total derivative) by set- 
ting lo\ = uj 2 = 0, which will be the case at the end of 
our calculation, hence the total linear momentum is con- 
served and the dynamics of the system cannot, therefore, 
change the position of the center of mass Q . Notice that 
the presence of the magnetic field makes the Lagrangian 
strictly two-dimensional. 



EFFECTIVE ACTION 

We proceed to the effect of an external transverse mag- 
netic field on the tunneling of the particles out of a 
metastable state by following the method of Caldeira and 
Leggett outlined in The partition function is given 
by 



Z 



dx\dx 2 K(x?\,x 2 , /3), 



where 



K(xi,X2,0) = J ' Vx[ J Vx 2 cx-p (-S E ) , 



and 



Sr. 



drCi 



(9) 
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where (3 = 1/T is the inverse temperature. The tunneling 
rate is proportional to cxp(— S E ) where the Euclidean 
classical action S E is determined by the bounce solution 
of the equation SSe = in which the periodic boundary 
condition x[(0) = x~i(/3) and 2^(0) = 22 (/3) are required. 
We shall set H = 1 throughout our calculation in this 
paper. Let's simplify the problem by taking mi = m 2 = 
m and uji = co 2 = oj' and now introduce the following 
change of variables 



q = xi- x 2 , 
x\ + x 2 



(12) 



where q is the position of particle 1 relative to particle 
2 and r is the position vector of the center of mass of 
particles 1 and 2. In the new coordinate system, the 
Lagrangian takes the form 



Ce = T^mqi 2 + V(\qi\) + ^rhuj' 2 q 2 

M 2 M /2 2 

+— n +yw ri +iqBxe ij r i q j , 



(13) 



where m = |m is the reduced mass and M — 2m is the 
total mass. We have dropped a total derivative in the 
Lagrangian which comes from the magnetic field term, it 
has no contribution to the classical equations of motion. 
The density matrix becomes 

rq(P)=<ib rr(P)=f 
K(q, r, (3) = Vq Vfexp (-S E ) ■ (14) 

Jq(0)=q Jf(0)=ro 

Exploiting the periodic boundary conditions on (f and r, 
one can expand them in terms of Fourier series [J] : 



n— 00 



n— — 00 
n— 00 



(15) 



where <f_„ = , r_„ = r* n and u> n = —u)- n = 2im/fi is 
the bosonic matsubara frequency. The classical equation 
of motion for r is 



Mn + iqB <,,(/', - Mu/Vj = 0. 
Fourier transforming (|16p we obtain 

qB±L0 n €ijqj n 



(16) 



(17) 



M{wl + lo' 2 ) ' 
For any path, we write 

n(r) = n(r) + w (t), (18) 
so yi(0) = Vi{P) = 0. The integral over fi(r) becomes 



: (j dT \^r 2 + ^Mu' 2 r 2 +iqB^ l] f l q J 
["drf^yf + hl^y 2 



(19) 



The linear terms in yi vanishes by virtue of the classical 
equation of motion (|16[) . Fourier transforming (|19p and 
using (|T7|) we obtain 



^ = 7» £ ~M(u 2 +uj> 2 )\y v 



n— — 00 
n—oo 



3 1 



{qBj_ojr, 



(20) 



p ^ 2M{loI+u' 2 ) 

n— — oo 



I Qin I 1 
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using the fact that Vri = T>n n — T>yi n , the Gaussian 
integral over is easily done, we finally obtain an ef- 
fective action given by 



jeff 







where 



2/3 



A n = - I rhuj 



^ n—oc 

7q ^ ] •^nlftnl i 



(21) 



{qBj_uj ri 



M(cj2 + u/ 2 ) 

The effective action can be written equivalently as 



(22) 



S«"=j^dr{|* 2 + V(|g«|) 



- I dr / dr^(r-r)gi(r)gi(r), 
^ Jo Jo 



(23) 



where 



^ n—oo 

Ar-r') = - E A » 



>(t-t') 



(24) 



The first terms in (|22[) is independent of w n and thus give 
an unnecessary delta function contribution to (|24[) which 
does not have any contribution to the effective action and 
hence can be neglected. This allows one to write ([24]) as 



/ ~ D \ 2 1 ™-°° , ,2 »W„(T— T' 



V M ; -(25) 



Using the fact that 



/ dT^(r) = 0, 
Jo 



(26) 



the equivalent form of the Caldeira Leggett effective ac- 
tion [l( one can obtain from (|2"3")l is 



S 



eff 



+ i f dr ! dr'A(r - r') [q t (r) - qi (r')} 2 



(27) 



Further simplification of (|25|) yields 



A(f) = M 



( qBj 



2 n— oo r 



\ M J P 
qB ± 



\ t 



M 



n— — oo 

2 n— oo 

a E 



(28) 



M J ft 



K +W '2) 



This allows one to use the summation formula 

= coBhju^/2-Jrj)] 
lTJ 2 V ill / sinh(/3w'/2) 

Exponentiating the hyperbolic functions we have 



(29) 



.A(r) 



Mw' ( qB x 



2 \ M 



(30) 



where f = r — r' and n(w) is a single-particle bose 
distribution function given by 



n(w') = 



1 



(31) 



We follow the method in [l| and periodically continue the 
path qi(r) outside the range < r < /3 by the prescrip- 
tion qi(r + 0) = qi(r), then (|3U|) becomes 



•4>(r) = ^ 



A/ 



(32) 



which is the zero temperature limit of (|30p . We can fur- 
ther simplify (|30|) by taking the limit w' — !> 0, this gives 



A 







M I qB ± 



M 



(33) 



The effective action becomes 



,5' 



eff 



dT{-e+v{\ qi \) 



M fqB d 



4/3 V M 



2 r P rP 

dr / dr' fe(r) - %(r')]' 
o Jo 



(34) 



It was shown by Ping [6j for one particle in three dimen- 
sions with a metastable potential only in the z-direction 
and a weak harmonic oscillator potential in x and y 
directions under the influence of an external magnetic 
field B = (Bn,0, B±) whose vector potential is given by 
A = (0, B±x — Buz, 0) that the effective action maps onto 
a one dimensional version of (|2"T|) and for Bu = which 
corresponds to Landau gauge, there is no dependence of 
effective action (tunneling rate) on B±. However, for the 
case of two charged coupled particles with equal and op- 
posite charges in the plane in the presence of an external 
potential and a traverse external magnetic field deter- 
mined by a symmetric gauge, and a metastable potential 
given by the interaction potential, the effective action is 
strictly two dimensional (i = 1,2) unlike the one dimen- 
sional version obtained in [lj and Q with the inclusion of 
an environmental coupling in the absence of an external 
magnetic field. 
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LEGGETT'S PRESCRIPTION AT ZERO 
TEMPERATURE 

In order to obtain the effective action at zero temper- 
ature (/3 — > oo), we appiy the Leggett prescription 
which tell us that if the Fourier transform of the real 
time classical equation of motion is of the form 



(^) 



(35) 



then the formula for the tunneling rate can be obtained 
from the effective action 

1 f°° 1 

S%" = ^ / du>-K(-i\u;\)\q(u>)\ 2 + S v (q(u>)), (36) 



where 



S„(<ZM) = / drV(q(T)). 



(37) 



At zero temperature, the Fourier series in real time takes 
the form 



(38) 



1 f°° 



Now the real time equations of motion from (|13j) are 

M(fi + oj' 2 ri ) - qB ± e ijqj = 0. (39) 



dV 
dq t ' 



(40) 



Next, we Fourier transform (|3"9"| and (|4"0)) and solve for 
qi(u>). The result is 



where 



K(cj) 



~ 2 - !'2 

-muj + muj 



M(w 2 - cu /2 )' 
Substituting (|4"2")l into (|3"6")l , we obtain 
1 Z" 00 1 

(gB^) 2 



(41) 



(42) 



(43) 



\qi(u))\ +S v (qi(uj)). 



M(w 2 W 2 ) 



In other to see how this result follows from plj) and ([22]) 
as /3 — )■ oo , we first Fourier transform the kinetic term of 
pip and replace u n by tu and the summation over n by 
integration over lo with a normalization factor of 1/2-k. 
Therefore, we see that the results are consistent. 



DISSIPATIVE ENVIRONMENT 

In this section we shall consider the coupling of the 
Lagrangian in to a thermal harmonic oscillators in 
two dimensions, therefore its classical dynamics will be 
dissipative. The Euclidean Lagrangian we shall consider 
is given by 



C E = —\x 1 \ +— \x 2 \ 



X\ X X\ — X 2 X X 2 



+ -m-i_ulx\ + \^m 2 u\x\ + V{\x! - x 2 \) / 44 j 



N 



a=l 



x 2 a +uj 2 a {x a ~ xY 



where x = {x\,x 2 ). 

In the absence of a transverse magnetic field and 
weak harmonic oscillator potentials, the Lagrangian cor- 
responds to a two dimensional version of the one con- 
sidered in Q except for the coupling in x instead of 
x 2 . The problem of one particle in three dimensions 
with a metastable potential only in the z-direction and 
a weak harmonic oscillator potential in x and y direc- 
tions under the influence of an external magnetic field 
B = (Sii, 0, B±) whose vector potential (Landua gauge) 
is given by A = (0, B±x — B\\z, 0), coupled to an environ- 
mental harmonic oscillators has been considered in 0|. It 
was shown that the problem can be mapped onto a one 
dimensional one. In (|44[) we have assumed that only par- 
ticle 2 is coupled to a large environmental harmonic os- 
cillators. Notice that the Lagrangian is still translational 
invariant up to a total derivative when u>i = uj 2 = 0. In 
order to obtain the effective action, we expand x a and x 
in a Fourier series as 



n— — oo 
n— oo 



*M = \ t 



(45) 



and perform the Gaussian integration over x an . This 
gives 

K{xi , af 2 , /?) = J Vxi (t) J Vx 2 (t) exp {-S E ) , (46) 
where 



Se = I dr 
h 



X\ x x\ — x 2 x x 2 



^(muj' 2 + ^2 m aU 2 a )x 2 + V(\x~i - x 2 \ 

a 

1 m a uj^ _ 2 

^ 2{ujl+uj2 fr t \ ■ 



(47) 
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We have set m\ — rri2 — m and uj\ — ui2 — w 1 to arrive 
at this result. Next, we rewrite (|47|) in terms of q and r 
using (fT2"j) , Fourier transform using (fT5|) and use the fact 
that the Fourier coefficient \x n \ 2 = |2?i n | 2 + |a?2n| 2 and 



(48) 
(49) 



Then the Gaussian integration over r ra can be easily done, 
and we obtain 



K(q,r;P) = J Vq(r) cxp (-S 



e J f ), 



(50) 



where 



S 



l ff = j Q dr^ 2 + V(\q\) 

* Jo Jo 
1 

+ - dr dr'B{T-T')q(r)q{T'), 
1 Jo Jo 



and 



n 



(52) 
(53) 



In this case the A n and B n are given by 



J\..,-\ 



(qB±u n f 



M( W 2 +W '2) + 2A„(/3) W 2' 



and 



B n = -A»C8). 



up) = E 



(54) 
(55) 



(56) 



We have dropped a constant term independent of uj n 
since it have no contribution to ()51|) . 

Now let's consider the limit of weak dissipation, that 
is A„(/3) << M(uj 2 + lu' 2 ), we have 



(qB±uj n ) 



M(u 2 n + u> 2 ) + 2u } 2 \ n (f3) M( w 2 W 2 ) ' 



(57) 



We then substitute (|5"7| into ([5"4"|) and use the definition 
of A n (/3) and the summation formula Q| 



1_ x - m a uj 2 a _ m a uj a cosh [uJJ3j2_ - f )] 



2 sinh(w Q /3/2) 



(58) 



we obtain in the limit ui 1 — > 0: 



+ - [°° dr' I dra(r - r') [q[r) - q{ T ')] 2 
4 J -oo Jo 



where 



y(t-t') = ^m Q ^cxp(-w a |r-T'|) (go) 



and uj c = qB± /M is the cyclotron frequency. 



OHMIC DISSIPATION 

In line with we shall assume that the effect of the 
oscillators on the motion of the particles result in the 
force of friction r)x. This requires that the spectral den- 
sity should be defined as 



(61) 



All the information concerning the effect of the environ- 
ment on the dynamics of the particles is contained in 
J (a)). We shall consider the case of ohmic dissipation in 
which J(uj) is given by 

J(u)) = r]ui, (62) 

up to a frequency cutoff uj c . Using (|6ip we have 



(63) 



7T In U) UJ 2 +UJ 2 



Substituting (|62[) into (|63p and performing the integra- 
tion, we obtain 



X n = r]/\uj n \. 



(64) 



Substituting this expression into the expressions for A n 
and B n and setting oj' = we have 



je// 







+ 7 f dT I dT ' A( - T - r ') W) - «V)] 2 (65) 
4 Jo Jo 

1 rP rP , 
+ -J dr J dr'B{r - r') [q(r) - q(r')f , 



where 



A( T - T >) - -1 V (l B ^n) 2 iUn (r-r') (m) 



(i 



Let us now work out the effective action at T = for the 
case in which the motion of the two identical particles are 
damped by Ohmic friction with constant rj respectively. 
We shall do this by applying the Leggett's prescription 
outline above. The real time classical equations of motion 
are 



where 



qB A 
qBj 



dV 
dx]- 
dV 
dx 2 ' 



(68) 
(69) 



Substituting x\ and x 2 in terms of and qi using the 
inverse transform of (|12[) . we obtain, after adding and 
subtracting the resulting equations 



2mf j - qB^eijijj 



1 



-mq\ - ql! «,,/', 



1 



o, 



:V9i 



dV 

■2 " ' L ."" 

Fourier transforming (|70|) and solving for qi(to) 



dV 

K(ui)qi(ui) = - — (uj), 
oqi 

where K(ui) in this case is given by 

~ 2 , {qB^ f 



(70) 

we obtain 
(71) 



K(lu) 



irjuj. 



-Mlu 2 + 2ir)U) 
The effective action at T = is then given by 



(72) 



S 



eff 
E 



1 

" 2^ 
2 



duj- 



|ft(w) 



fhu! 2 



(qB±u) 



Muj 2 
S»(?<(w)). 



2r)\i 



(73) 



We also noticed that the effective action found in [Icj for 
the case of no pinning and finite dissipation at T = does 
not give the same expression using Leggett's prescription 
at T = for the case of Ohmic dissipation. 

Now consider the limit B± — 0, in this limit, the ex- 
plicit expression of (|57|) can be easily obtained. The ef- 
fective action becomes 



s e J f 



13 





8tt 



(It 





fhq 2 



dr' 



dr- 



[q{T)-q{r>)Y 



(74) 



This effective action is similar to a two dimensional ver- 
sion of the mi = ?Ti2 limit of Eq. (|44| in Ref.Q. Now, 
consider the limit of strong dissipation Mu^ << 2rj\uj n \ ) 
in this limit, the explicit expressions in (|66|) and (|67p can 
be summed. The effective action becomes 



S 



eff 



dr 



in 

Veff 

8tt 



-rhq 2 



V(\q\) 



dr' 



[q(T)~q(T>)} 2 
o - (/3/7r) 2 sin 2 (^//3)' 



(75) 



dr- 



Veff = 



(qB±f 



(76) 



CONCLUSIONS 



We have studied the macroscopic quantum tunneling 
of two coupled particles of identical masses but opposite 
charges in the plane, in the presence of a constant trans- 
verse external magnetic field and an external potential 
whose interaction potential allows for a metastable state. 
We showed that the effect of the magnetic is to suppress 
the tunneling of the particles out of a metastable state 
and also our effective action remains two dimensional un- 
like a one dimension version obtained in We further 
coupled the system to a thermal bath of harmonic oscil- 
lator and showed that in the limit of weak dissipation, 
there is an effect of the magnetic field to the effective 
action. In the limit of zero magnetic field, we repro- 
duced a two dimensional version of the results obtained 
in Q for mi = mi, which also coincides with the results 
of Caldeira and Leggett [lj in two dimensions in terms 
of the reduced mass. Our results are as a consequence 
of the conservation of total linear momentum as was first 
shown in Q since the magnetic field term is translational 
invariant up to a total derivative. The results obtained 
here can also be applied to the metastable states of the 
molecules of fluid and solids in a transverse magnetic field 
and coupled superconductor vortex tunneling in two di- 
mensions. 
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